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Abstract 

We performed molecular statics calculations of energy barriers for adatom moves 
in the vicinity of steps on Pt(lll) surface. We used the semi-empirical many-body 
Rosato-Guillope-Legrand potential and we systematically calculated barriers for 
descent of straight steps, steps with a kink and small islands as well as barriers for 
diffusion along the step edges. We confirmed that the lowest barrier for descent 
is for an exchange process near kink's or island's corner on a step with a {111} 
microfacet. Diffusion along a step with a {111} microfacet is faster than along a 
step with a {100} microfacet. We also calculated barriers for diffusion on several 
surfaces vicinal to Pt(lll). Our results are compared with previous calculations. 

1 Introduction 

The diffusion in Pt/Pt(lll) system has been recently subject of several studies both 
experimental [1-3] and theoretical [1,4-11] It was motivated by growth experiments, par- 
ticularly by observation of the so-called reentrant layer-by-layer growth |12| and also by 
growth-temperature-induced switch from a triangular growth shape bounded by steps 
with a {100} microfacet to one bounded by steps with a {111} microfacet | 13| . In order 



to explain these experiments on the microscopic level a detailed understanding of the 
diffusion on inhomogeneous Pt(lll) surfaces is needed. 

The self-diffusion energy barrier on the flat Pt(lll) surface has been determined accu- 
rately in two independent experiments. Bott et al. Q deduced from scanning-tunneling- 
microscopy (STM) measurements of the density of platinum islands at different tempera- 
tures an activation energy 0.26± 0.01 eV. This is in remarkable agreement with the recent 
field-ion-microscope (FIM) experiment [|3] giving the barrier 0.260± 0.003 eV. The agree- 
ment of theoretical calculations with the experimental value is rather unsatisfactory (see 
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Table |l|). Semi-empirical potentials systematically underestimate the experimental value 
and the existing ab-initio full potential (FP) calculations give in the contrary significantly 
higher energies. 

Energy barriers for adatom moves on an inhomogeneous surface are not easily ac- 
cessible by experiment and it is the aim of the theory to calculate energy barriers for 
elementary processes which are supposed to be relevant for explanation of growth ex- 
periments. Several calculations have been already performed. They confirmed that the 
layer-by- layer growth at low temperatures can be explained by decrease of step-edge bar- 
rier (called Ehrlich-Schwoebel barrier) for small and irregular islands || |fj, ||. As regards 
the shape of growing islands it has been analyzed by several theoretical groups J|, [7], || , 
but so far no mechanism has been agreed on. 
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Table 1: Results for self-diffusion barrier Eg (in eV) on flat Pt(lll) surface, EAM - 
embedded atom method, EMT - effective medium theory, FP - full potential calculations, 
for other abbreviations see text. 



In the situation when there is not quantitative agreement between theory and experi- 
ment already for diffusion on the flat surface, one can still look at trends and compare mu- 
tually values of energy barriers for various adatom moves obtained within one approach. 
However, different authors investigated various processes using different potentials and 
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the comparison is not always possible. Therefore, it is desirable to compare barriers for 
as many as possible processes using one approach. 

In this paper we use the semi-empirical many-body Rosato-Guillope-Legrand (RGL) 
potential [ |TS| , |T9"[ developed within the second moment approximation of tight-binding 
theory. This potential allows relatively fast calculations and, at the same time, provides 



reasonably reliable results pOfl . We show that with the RGL potential we can reproduce 
majority of previous results. In particular, we systematically study exchange processes 
for descent of steps and diffusion along straight step edges. 

2 Method 

The simulations were done with finite atomic slabs with a free surface on the top, two 
atomic layers fixed on the bottom, and periodic boundary conditions in the two directions 
parallel to the surface. The slab representing the substrate of (111) surface was 11 layers 
thick with 448 atoms per layer. For diffusion along channels on the vicinal surfaces (311), 
(211), (331), (221) and (322) we used the systems of approximately 5000 atoms that 
contained 19 to 44 layers, with 110 to 240 atoms per layer. 

The equations of motion were solved using a leap-frog algorithm with a time 
step of 0.5 • 10~ 14 s. We used the classical NVE ensemble and molecular-dynamics cooling 
method for molecular statics calculations of energy barriers. A conventional spherical cut- 
off with the cutoff radius of 6.4 A was used. We have applied the RGL potentials with 
the parameters obtained by a fit to experimental data taking into account interactions up 
to the fifth nearest neighbors |L9] . 



The energy barrier of a particular diffusion process was obtained by testing system- 
atically various paths of adatom moves, and the path with the lowest energy barrier 
was chosen to be the optimum path. The adatom diffusion barrier Ed is defined as 
Ed = E sa d — E m i n where E sa d and E m i n are the total energy of the system with the 
adatom at the saddle point and at the equilibrium adsorption site, respectively. The 
minimum energy paths for both the jump and exchange processes were calculated. The 
minimum energy path for jump diffusion was determined by allowing the migrating atom 
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to relax in a plane perpendicular to the path at each step. The rest of the atoms in the 
system was allowed to relax in all directions. The energy barrier of the exchange process 
was obtained by moving the surface atom, which was to be replaced, by an adatom with 
finite steps along the direction of exchange. This atom was allowed to relax in the plane 
perpendicular to the exchange direction at each step, whereas the other atoms, including 
the adatom, were allowed to relax in all directions. 

3 Results 

3.1 Descent to the lower terrace 

As the first step we calculated the energy barrier for diffusion on the flat Pt(lll) surface. 
We found the value 0.17 eV, which is comparable with the values obtained from the 
effective medium theory (see Table [I]). 

We studied descent of an adatom to the lower terrace from both types of steps on 
the (111) surface, i.e. , step A with {100} microfacet and step B with {111} microfacet 
(see Fig. 1). We performed calculations for several different geometries: straight steps, 
steps with a kink, island 3x3 and island 6x6. We considered both direct jumps of an 
adatom and exchange processes. For all considered geometries we systematically studied 
all possible adatom moves, and investigated the dependence of energy barriers on the 
local geometry. Our results are summarized and compared with previous calculations in 
Table |. 

Let us consider at first the case of straight steps. The energy barrier for direct jump 
from the upper to the lower terrace is 0.53 eV for step A and 0.52 eV for step B. The 
barrier for exchange at the step A (0.50 eV) is comparable with the barrier for jump, but 
at the step B, the barrier for exchange (0.40 eV) is significantly lower. We did not find 
a difference in binding energy between fee and hep site on the flat surface. However, we 
observed a difference between binding energies at fee and hep positions near the steps 
on the upper terrace. We found that near the step A the fee site is energetically favored 
by 0.02 eV in comparison with the hep site, but on step B the hep site is energetically 
favored by the same value in comparison with the fee site. 
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Table 2: Energy barriers (in eV) for descent of steps on Pt(lll). For the meaning of 
labels see Fig. 2. 
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Table 3: Summary of energy barriers for exchange on a kink on step B. 
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In geometry with a kink on a step one has to consider more processes. As an example 
of our approach we show in Table ^ the energy barriers for different exchange processes 
near a kink on step B. We calculated energy barriers for all possible exchange events and 
for each exchange process all moves of the replaced atom were considered. For example 
in the case of the exchange process presented in the first line of Table |3| (the adatom 
starts in the fee site labeled by 1 and pushes out the edge atom rfl) two possible different 
paths schematically shown by arrows in Fig. 4 were considered. Only the minimal barrier 
heights are presented in Table [| 

There is only small influence of kink geometry on energy barrier for direct jump (Tab. 
K). Contrary to Liu and Metiu |§ but in agreement with Villarba et al. |J and Jacobsen et 
al. ]7) we found that the diffusion barrier for exchange on steps with kinks is significantly 
lower than on straight steps. We found that the smallest energy barrier for exchange 
on the step with a kink is on the step A for exchange over the kink atom (exchange of 
adatom from the hep position 5 with the edge atom rf3 - see Fig. 2) and on the step B 
for exchange next to corner (from the hep position 2 to the edge atom position rfl). 

We continued the investigation of geometry effects on energy barriers for descent with 
the study of small islands (3x3 and 6x6 atoms) on Pt(lll) surface. We found no 
changes in the barrier heights for direct jumps. For exchange processes on small island 
further reduction of energy barriers has been observed. For 3x3 island minimal values 
were obtained for the exchange of the atom in the middle of the edge (0.26 eV for A-type 
edge and 0.13 eV for B-type edge). For the 36-atom island the minimal energy barriers 
(0.30 eV for A-type edge and 0.16 eV for B-type edge) were found for the exchange pro- 
cesses starting in the hep sites in the vicinity of the island's corner moving away the edge 
neighbors of the corner atoms. The exchange mechanism is more important for diffusion 
on small islands than for diffusion in the vicinity of longer atomic edge. 
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3.2 Diffusion along the step edges 

We also calculated barriers for diffusion along straight steps (Fig. 3). We found that 
diffusion along the step of type A is slower (the barrier is 0.55 eV) than along the step 
of type B (the barrier is 0.43 eV) - see Fig. 3. This is in agreement with several previous 
calculations (0.44 eV and 0.40 eV) g, (0.60 eV and 0.43 eV) g, (0.23 eV and 0.18 eV) 
0. But it is in contradiction with the assumption made in the recent kinetic Monte Carlo 
simulation || attempting to explain island shapes observed in homo-epitaxial growth of 
Pt(lll). The authors justified their assumption by comparison with the results of field- 
ion-microscope measurements on the (311) and (331) surfaces ]|T4"| , p2| . These surfaces 
consist of narrow terraces with step edges of type A and B, respectively. Experimental 



values of energy barriers are E 3 u = 0.60 ± 0.03 eV ||22|| , and E 3 u = 0.84 ± 0.1 eV 
TA\. Nevertheless, it may be misleading to use barriers obtained for the (311) and (331) 



surfaces instead of barriers along the step edges on the (111) surface. It is apparent (see 
Fig. 4) that when not only nearest neighbors, but also more distant neighborhood is 
considered, then the geometry is different and one can expect different contributions of 
adatom-surface interactions to the energy barriers. The relaxation of vicinal surfaces is 
also different. In principle one cannot exclude even qualitative difference as it seems to 
be the case for platinum and as it has been observed on iridium [p3[] . Unfortunately up 
to our best knowledge, there is no experimental result for diffusion along step edges on 
Pt(lll). 
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Table 4: Comparison of activation energy barriers (in eV) for diffusion on the (311) and 
(331) surfaces. 



Table [| shows results of previous calculations for diffusion along steps on the (311) 
and (331) surfaces. Three EAM calculations in Table |] differ by parameterization. We 
can see that all EAM calculations give a lower barrier on the (331) surface than on the 
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(311) surface which is in contradiction with the experiment. Only the calculation using 
the Morse potential gives qualitative agreement with the experiment. 

In our calculation using the RGL potential, we obtained the same barriers 0.50 eV 
for diffusion along steps on both (311) and (331) surfaces. We calculated also barriers 
for diffusion along steps on some other vicinal surfaces with wider terraces: (211), (221) 
and (332) shown in Fig. 4. All our results are summarized in Table || There is a clear 
tendency with the changing distance between steps. In the case of surface geometries 
with A-type steps the energy barrier is increasing with the step-step distance whereas for 
the B-type steps it is decreasing. We note that the same trend is valid also for the data 
reported in Ref. ||. 
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Table 5: Activation energy barriers (in eV) for diffusion along steps on different surfaces 
calculated with the RGL potential. 

4 Conclusion 

We systematically calculated barriers for descent of straight steps, steps with a kink and 
small islands, and we determined the processes with a minimal energy barrier. We have 
shown that with the many-body RGL potential we can obtain results comparable with 
those obtained from the effective medium theory. We found that the descent is easier at 
the step B than at the step A and that for irregular geometry the exchange mechanism of 
diffusion must be taken systematically into account. The lowest energy barrier for steps 
with kinks is significantly lower than for straight steps and the energy barrier for descent 
from small islands is even much lower. 

In the case of diffusion along step edges we obtained lower barriers for diffusion along 
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the step B than for diffusion along the step A. Calculations for vicinal surfaces demonstrate 
the importance of the step-step interaction on the energy barriers. We observed the trend 
that the diffusion along the step B on Pt(lll) surface is faster than along the step A. 
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Figure captions 



Fig. 1. Two types of step edges, A and B, for a large island. Solid line shows the 
diffusion path along the island edge. The atomic layers from the surface to the bulk are 
large filled circles, large open circles, small open circles, and very small open circles. 

Fig. 2. Different exchange processes near a kink site on the step B on Pt(lll) sur- 
face. The edge atoms undergoing investigated exchange diffusion (rfl,...,rf4) and starting 
positions of an adatom (1,...,9) are shown. 

Fig. 3. The dependence of adatom energy for the path along the edge of a large is- 
land shown in Fig. 1. 

Fig. 4. Surface normal views for Pt surfaces where self-diffusion along the A-step type 
edge, i.e. (311) and (211), and along the B-step type edge, i.e. (331), (221) and (332), 
has been calculated. Solid line shows the diffusion path along the edge for each surface. 
The smaller filled circles are, the deeper the atomic layers lie in the bulk. 
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